Strong Dynamical Heterogeneities in the Violation of the Fluctuation-Dissipation Theorem in Spin 

Glasses 
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We analyze numerically the violation of the fluctuation-dissipation theorem (FDT) in the ±7 Edwards- 
Anderson (EA) spin glass model. Using single spin probability densities we reveal the presence of strong 
dynamical heterogeneities, which correlate with ground state information. The physical interpretation of the 
results shows that the spins in the EA model can be divided in two sets. In 3D, one set forms a compact 
structure which presents a coarsening-like behavior with its characteristic violation of the FDT, while the other 
asymptotically follows the FDT. Finally, we compare the dynamical behavior observed in 3D with 2D. 

PACS numbers: 75.10.Nr, 75.40.Gb, 75.40.Mg 
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The study of glassy behavior, characterized by out-of- 
equilibrium dynamics that present very long relaxation times, 
involves a vast range of different systems such as spin glasses 
(r], structural glasses, colloidal and polymeric systems, and 
granular systems to name just a few [2]. These systems 
present out-of-equilibrium properties such as aging and vi- 
olation of the fluctuation-dissipation theorem (FDT) Jsl 01, 
whose characterization can be used as signatures of typical 
dynamical behaviors. For example, according to their devia- 
tion from the FDT they can be classified into three different 
groups: coarsening, structural glass and spin glass systems 
||31]. These cases have been theoretically described in terms 
of replica-symmetry breaking (RSB): in coarsening systems 
replica symmetry is unbroken, structural glasses correspond 
to one-step RSB and spin glasses to full RSB |30. In the case 
of short-range spin glasses, most results on the violation of 
the FDT have been qualitatively and quantitatively described 
in the framework of mean field theory (full RSB) yilsl]. Also 
direct experimental evidence on the violation of the FDT in an 
insulating spin glass |6] have been fitted by mean field theory. 

In this work we show that, in contrast to these results, the 
violation of the FDT in the 3D ±7 Edwards-Anderson (EA) 
model is the result of two components with completely dif- 
ferent behaviors: one that tends to satisfy the FDT relation, 
and another which presents a violation of this relation similar 
to coarsening systems. The behavior of the latter component 
seems to be compatible with the droplet picture scenario ■ 

We reveal the presence of two components in the viola- 
tion of the FDT by a careful analysis of dynamical hetero- 
geneities. In principle, two different approaches may be used 
to study the development of dynamical heterogeneities. On 
one hand, it is possible to use space or time coarse-grained 
protocols |8, ^. On the other hand, it is possible to use sin- 
gle spin observables looking for a direct observation of the 
local heterogeneities. In any of these approaches, if one is 
interested in making disorder averages and chooses to iden- 
tify the spins by their position in the lattice, trivial results are 
obtained, since the difference between spins are washed out 
ifioll . Here we take a local approach, but in contrast with pre- 
vious works, we properly include disorder averages by using 



a constrained structure of the ground state called backbone 
(see below), which allows for much insight in the analysis of 
dynamical heterogeneities Ol ill . 

We begin our analysis by using single spin probability den- 
sities to show that dynamical heterogeneities are present in 
the violation of FDT in the 3D ±J EA model. These hetero- 
geneities present a non-trivial complex structure with bimodal 
distributions 111. 12.,, 13.1 , and we refer them as strong, in con- 
trast to weak dynamical heterogeneities which present uni- 
modal distributions with an elongated shape. Then, as done 
in 2D in 1 11], we use information from the ground state topol- 
ogy to establish a quantitative relation between spatial and dy- 
namical heterogeneities. This analysis allows us to present a 
new physical interpretation of the violation of FDT. Finally 
we compare the violation of FDT in 2D and 3D, which re- 
veals fundamental differences. 

We consider the EA model for spin glasses UJ, defined on a 
D-dimensional lattice of linear size L with periodic boundary 
conditions. The Hamiltonian of the model is 



H 



(iJ) 



(1) 



where cr, = ±1 is the spin variable and indicates a sum 
over the 2D nearest neighbors. The coupling constants Jjj = 
±7 are random variables chosen from a symmetric bimodal 
distribution. The time evolution of the model is governed by 
a standard Glauber dynamics with sequential random updates 
using a continuous time Monte Carlo algorithm fT?]. 

In order to investigate the violation of the FDT we use a 
two-times protocol which emphasizes the out-of-equilibrium 
character of the dynamics. The initial condition corresponds 
to a quench from T — °o to the temperature of interest T < 

= 1 . 12 lITsil at f = 0, where each spin takes a random value 
(7,- = ±1. From this initial condition different two-times quan- 
tities are analyzed, which depend on both the waiting time t„, 
when the measurement begins, and a given time f > [!3i,4]- 

We focus on the single spin two-times correlation function, 
defined as % = {Ci{t,t^^)) — {ai{t)ai{tw)) , with t > t^, and 
where (. . .) is an average over thermal histories, that is over 
different initial conditions and realization of the thermal noise. 
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In terms of the single spin correlation, the global two-times 
correlation function can be expressed as 



1=1 



(2) 



where N = LP and [...]„,. indicates average over different re- 
alizations of bond disorder (samples). In the same way the 
global integrated response function is defined in term of the 
single spin integrated response function as 



(3) 



where Xi = {'^i{t)''>ign{hi{t„)) / h) and, as usuallv\ a random 
field of intensity h is switched on at time t„ IB |4|] . All the 
results presented correspond to T = 0.8, h = Q. \ and 10^ ther- 
mal histories for each sample. 

Using the global quantities defined in Eqs. (|2]l and 
the out-of-equilibrium fluctuation-dissipation relation can be 
written as 



Tx{t,t„)^X{C)[\-C{t,u)] 



(4) 



where X is the fluctuation-dissipation ratio (FDR) d. A 
useful representation of Eq. ^ is the parametric plot of 7% 
vs. C iS H. When the FDT holds the FDR is X = 1 and 
the parametric plot shows a linear relation with unitary slope. 
In an out-of-equilibrium situation the FDT does not longer 
hold and two regimes are observed: for t / <C 1 the system 
shows quasi-equilibrium with X ~\, while for f/f,v 3> 1 a vi- 
olation of the FDT is observed with X <\. The behavior of 
X{C) for f/f,v ^ 1 allows for a simple classification of out-of- 
equilibrium systems into three main categories: (/) the value 
X = is related to coarsening systems, (//) a constant X <\ 
value is associated with structural glasses, and (///) a decreas- 
ing monotonic < X{C) < 1 function is associated to spin 
glasses. We will later use this simple scheme to physically 
reinterpret the out-of-equilibrium dynamics of spin glasses. In 
particular, the continuous line in Fig. [1] shows the parametric 
plot for the 3D EA model with linear size L = 20 for = 10. 
The results correspond to averages over 50 different samples. 
After the departure from the quasi-equilibrium regime a con- 
tinuous variation of the FDR is observed 1 12i [l7ll . 

In order to study how the heterogeneities arise in the viola- 
tion of the FDT, we measure the joint probability distribution 
(JPD) of the single spin quantities, p['i^(f,fvv),r;^,(f,f„,)]. In 
Fig. [T|we show the projection of p in the {Tx , C) plane for the 
same 50 samples used above. The sequence shows the evolu- 
tion of the JPD for increasing At — t — t^, as a map plot. The 
parametric plot (solid curve) of the global quantities is indi- 
cated in each figure as a guide to the eye, and for each value 
of At the mean value of the distribution is indicated with a 
white circle. 

Note that already for very short Af the distribution presents 
an elongated shape indicating an heterogeneous behavior As 
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FIG. 1: Projection of the JPD of single spin correlations and inte- 
grated responses, p[%{t,t„),TXi{t,tw)], for L = 20 and f,„ = 10. The 
parametric plot of the global quantities, Tx{C), is reproduced in all 
the figures as a continuous line. The white dot indicates the mean 
value of the distribution for each At. 



At grows the distribution elongates further and the presence of 
two peaks can be clearly observed. At this point, it is clear that 
the mean value of the bimodal distribution shown in Fig.lTJd) 
cannot be used as a representative quantity to characterize the 
parametric plot. Since for the development of the two peaks 
it is necessary to reach the longest At as possible, we chose 
to work with a low t^. = 10 value. For longer t^ values we 
observed the same trend to a bimodal JPD for long time dif- 
ferences. 

In the search for the origin of these strong dynamical het- 
erogeneities we focus our analysis on a constrained structure, 
which is usually called a backbone. The presence of a back- 
bone in models defined on randomly connected graphs leads 
to time-scale separation and heterogeneous dynamics at suf- 
ficiently low temperatures jlSll . A backbone is also present 
in the EA model 23], which is formed by those bonds 
which are always satisfied or always frustrated in all ground- 
state configurations. This structure can be used to divide the 
spins in the system in two sets: solidary spins, which are the 
spins in the backbone and thus maintain their relative orien- 
tation in all the ground-state configurations, and non-solidary 
spins which are simply all the spins that are not in the back- 
bone 111 in . Further insight in the static properties of the 3D 
and 2D EA model has been recently addressed by Roma et 
al. 12 ll I22II . They found that the fraction of solidary spins 
in 3D(2D) is 0.76(0^7) of die total. In contrast with what 



is observed in 2D is El, the largest component of the sol- 
idary spins in 3D does not fragment as the system size grows. 
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FIG. 2: (a) Parametric plot for L = 8 (continuous line) and L = 20 
(open diamond) for five different waiting times. The arrow indicates 
increasing t„. (b) Projection of p and parametric plot of the global 
quantities Tx{C) for Z, = 8 (circle). We also present the separation 
into sets A (square) and B (triangle), (c) and (d) show p and its cor- 
responding parametric plot constrained to sets A and B respectively. 
In (b), (c) and (d), the white symbols indicate the mean value of the 
JPD distribution for = 10 and A/ = 10"*. 



involving 60% of all the spins. Furthermore, the largest com- 
ponent in 3D forms a compact percolating cluster Ii22ll . We 
use this information to analyze the strong dynamical hetero- 
geneities observed in Fig.[Tl and divide the system in two sets: 
the spins belonging to the percolating cluster of solidary spins 
(set A) and the remaining ones (set B). 

Since we are interested in performing averages over dis- 
order using this division, a considerable number of different 
realizations, where ground-state information is used, need to 
be taken into account. In order to calculate the backbone for 
each sample, we used an improvement of the Rigid Lattice 
Searching Algorithm as in ll22ll . This requires a large compu- 
tational effort, since in order to obtain the exact backbone it 
is necessary to reach the ground state 0{N) times. This sets a 
limit to the largest size we can consider to = 8^ ll22ll . 

To establish finite size effects we compare the results ob- 
tained with L = 8 (over 300 samples) and L = 20 (over 50 
samples). In Fig. |2ja) we show the parametric plot for L = 8 
(continuous line) and L = 20 (open diamond) for five differ- 
ent waiting times. For both system sizes the results overlap, 
which shows that, up to the times considered, no significant 
deviations due to finite size effects are present. Note also the 
shape of the JPD for L = 8 (Fig. I^b)) and L = 20 (Fig. [lid)), 
which shows that for both system sizes the qualitative results 
obtained are similar. 

Figure |2|b) also shows the behavior of the global paramet- 



ric plot (circle) and its division into sets A (square) and B 
(triangle). The JPD constrained to sets A and B are presented 
in Figs.|2lc) and|2ld) respectively. Note that the mean value 
of the JPD for the set A (B) is very close to the lower (upper) 
peak of the total JPD. In both cases the distributions are very 
broad, however the contribution of each set to the total distri- 
bution is clearly divided. Spins in set A (B) mainly contribute 
to the lower (upper) peak. 

We now present and support a simple physical interpreta- 
tion to understand the different behavior of each set of spins. 
Let us first highlight a significant topological property of set 
A. We observe that in average approximately 10% of the 
bonds in this set are frustrated in the ground state. By using a 
gauge transformation for a particular sample, which leaves the 
Hamiltonian invariant [23], one can change the bonds in set A 
so that all its spins are aligned in the ground state. As a conse- 
quence, this set is mapped into a ferromagnet with a low con- 
centration (« 10%) of anti-ferromagnetic bonds. Note that, 
when the anti-ferromagnetic bond concentration x is changed 
in the 3D EA model, it is well established that a ferromagentic 
order is present below x = 0.22 [24]. Then, since set A is be- 
low this limit (x w 0.1), we expect a ferromagnetic-like order 
to persist in this region up to a finite temperature. This result 
allows us to establish an unexpected property of the 3D ±7 
EA model: restricted to set A, an out-of-equilibrium behav- 
ior in accordance with that observed in a coarsening system 
should be found. 

In order to analyze the violation of the FDT, we present 
in Fig. [3] (a) the parametric plot for each set for ten different 
t„. We find that the lower curves, which correspond to set A, 
present a clear saturation to a curve away from the FDT line 
for increasing ?„., which is not observed when the system is 
considered as a whole (Fig.[3](b)). This behavior reflects the 
particular topological characteristics of set A. We expect that 
for larger system sizes and long times the curves will follow 
an out-of-equilibrium behavior similar to the violation of the 
FDT in coarsening systems [3]. 

In contrast to the saturation behavior of set A, the upper 
curves in Fig.[3](a) show that set B tends to the FDT limit for 
increasing f,,,. This surprising result is completely unexpected, 
since the out-of-equilibrium dynamics of the system is taking 
place below the critical temperature Tc, and a violation of the 
FDT should be observed when starting from random initial 
conditions (T ^ Tc) [2, J,]. Notice that all the information on 
the underlying processes, which we observe using the division 
into sets A and B, is lost when one considers the system as a 
whole (Fig.[3](b)). 

In order to stress the fundamental differences that can be 
pointed out in the intuitive frame proposed, we compare the 
out-of-equilibrium dynamics of the EA model in 2D and 
3D. These systems present very similar behavior at low tem- 
pera tures, allowing for comparisons between their dynamics 
till m 111. However, the fact that Tc = in 2D and Tc > 
in 3D, recalls that the physical mechanisms present in the dy- 
namics are clearly different. Notice also that since there is no 
percolating cluster of solidary spins in 2D [2Q], the system can 
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suits suggest that in short-range spin glasses the droplet pic- 
ture could be valid if restricted to a finite fraction of the sys- 
tem. However, as we have shown in the 3D EA model, to 
understand the behavior of the system as a whole it is also 
neccesary to take into account the existence of a finite frac- 
tion of the system which presents a paramagnetic behavior 
We believe that the inclusion of these elements in a new the- 
oretical description can contribute to the understanding of the 
physical nature of spin glasses. In particular it can contribute 
to the study of the relation between the percolation problem 
and the spin-glass transition 1, 261 . and also to hot topics such 
as temperature and disorder chaos ll27ll . 
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FIG. 3: Parametric plot for: (a) sets A and B and (b) for the system 
as a whole. The arrow indicates increasing t„ for the ten different 
waiting times. 

only be divided in solidary and non-solidary spins ifTlll . Using 
this division for the 2D EA model we found that, as expected, 
the evolution of the parametric plot for both sets tends to the 
FDT limit for increasing t^, in contrast to the results obtained 
in 3D. 

Summarizing, we observe strong dynamical heterogeneities 
in the JPD of single spin correlation and integrated response 
functions. We have gone beyond previous works on single 
spin quantities, and have been able to show that these hetero- 
geneities persist even when averages over disorder are taken 
into account. Using information from the ground state topol- 
ogy we are able to establish a quantitative relation between 
spatial and dynamical heterogeneities. In particular, we ob- 
serve that in the 3D EA model there is a fraction of the sys- 
tem which is capable of sustaining a ferromagnetic-like order. 
The presence of this structure is reflected in the parametric 
plot of Tx vs. C. On the other hand, the other fraction of 
the system presents a dynamical behavior that asymptotically 
tends to satisfy the FDT. As a consequence, the global out- 
of-equilibrium dynamics of the system can be regarded as a 
combination of two different components. A comparison of 
the violation of the FDT in 3D with 2D reveals that the under- 
lying physical mechanisms present in the dynamics are clearly 
different. 

We have shown that the continuous violation of the FDT, 
considered a signature of full RSB, can be regarded as the 
combination of two different components. This particular re- 
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